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Abstract 

It has recently been shown that the graviton can consistently gain a constant mass without 
introducing the Boulware-Deser ghost. We propose a gravity model where the graviton mass is 
set by a scalar field and prove that this model is free of the Boulware-Deser ghost by analyzing 
its constraint system and showing that two constraints arise. We also initiate the study of the 
model's cosmic background evolution and tentatively discuss possible cosmological implications of 
this model. In particular, we consider a simple scenario where the scalar field setting the graviton 
mass is identified with the inflaton and the graviton mass evolves from a high to a low energy scale, 
giving rise to the current cosmic acceleration. 



* huangqg@itp .ac.cn 
Vspiao@gucas.ac.cn 
' zhou.sy234@gmail.com 



1 



I. INTRODUCTION 



Although known as totally redundant degrees of freedom, gauge symmetries and general 
covariance are important for the modern understanding of particle physics and gravity re- 
spectively. While gauge symmetries can be straightforwardly broken to promote a massless 
gauge particle to a massive one, the equivalent for general covariance is more difficult to 
achieve Phenomenologically, the Standard Model gauge symmetries are spontaneously 
broken at low energies, where we are living, and gauge particles gain masses from the con- 
densate of the Higgs boson. It is interesting to ask whether a similar situation applies to 
the graviton, with the graviton mass smaller than the bound imposed by experiments? The 
recently observed cosmic acceleration may add another motive for considering this scenario 
(or modified gravity theories in general [2|), where the dark energy scale is introduced via 
the graviton mass. 

The linear massive spin-2 theory was discovered as early as Fierz-Pauli jjjj]. However, the 
gravitational force described by this theory differs significantly from that of linearized Gen- 
eral Relativity, for example, in the solar system. One can compute the scattering amplitude 
of Fierz-Pauli theory between two matter sources and find that when the mass goes to zero 
the amplitude generally does not reduce to the linearized General Relativity case. This is 
known as the vDVZ discontinuity |H, [EJ and seems to suggest that by simply performing 
experiments in localized regions one can infer that the graviton mass is mathematically 
zero, contradicting the expectation that forces mediated by bosons decay a la Yukawa. As 
pointed out by Vainshtein 0, in nonlinearized Fierz-Pauli models, the gravitational force 
should mimic the Einstein case in conventional environments for a small mass: A key ob- 
servation was that the introduction of the graviton mass has introduced a second nonlinear 
scale different from the one associated with the Planck scale. For example, around a spher- 
ical compact source, General Relativity is recovered within the nonlinear length scale called 
the Vainshtein radius, while linear Fierz-Pauli theory applies outside the Vainshtein radius. 
The Vainshtein radius is very large for a small graviton mass and goes to infinity when the 
mass goes to zero. However, nonlinear Fierz-Pauli models are not free from problems. The 
most pressing one is that there is a ghost originally identified by Boulware-Deser (BD) [tJ. 
Whilst a spin-2 particle should only have 5 propagating modes, nonlinear Fierz-Pauli gravity 
has 6 modes: The phase space of a generic 4 dimensional metric theory is 12 dimensional, 
as there are, after a Arnowitt-Deser-Misner (ADM) decomposition, 6 apparent dynamical 
components in the metric (the spatial ones); In nonlinear Fierz-Pauli gravity, there are no 
constraints at all, comparing to 8 constraints in the General Relativity case and 2 constraints 
in the linear Fierz-Pauli case. The sixth mode in nonlinear Fierz-Pauli models is responsible 
for the BD ghost @. 

Although broken in massive gravity, general covariance can be re-introduced by the 
Stiickelberg trick, i.e., introducing extra scalar fields patterned after general coordinate 
transformation [9|. In this way, we can easily decompose all the dynamical modes into a 
scalar, vector and tensor, and view the theory in the effective field theory's language, where 
the peculiarities of massive gravity such as the vDVZ discontinuity, Vainshtein mechanism 
and BD ghost can be clearly understood [of. In particular, most clearly seen in the decou- 
pling limit, the scalar sector contains higher derivative interactions that are strongly coupled 
at A 5 = (m^Mp) 1 / 5 (m g is the graviton mass.), which is rather small for a small mass. The 
scalar's equation of motion contains fourth order derivatives, so by a Cauchy surface dimen- 
sion counting, the scalar actually propagates two modes and, according to Ostrogradski's 
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theorem, one of the scalar modes is a ghost. 

However, as recently discovered, all the A5 
and the strong couplin g sc ale can be raised to A 3 = (m; 
Tolley massive gravity eliminating the BD ghost at least in the decoupling limit. 

This construction is closely connected to the galileon construction, which had been studied 



interactions can be systematically removed 
jMp) 1//3 in de Rham-Gabadadze- 



in a different class of modified gravity models (see e.g. [13Hl5|). In fact, the sixth BD mode is 



eliminated in the full dRGT model. This has been shown by analyzing the constraint system 



of the theory in the Hamiltonian formulation [10|, |12. 16- Unlike nonlinear Fierz-Pauli 
massive gravity, the dRGT model maintains the Hamiltonian constraint and a secon dary 
constraint is generated by the time conservation of the Hamiltonian constraint [l6|, 17 



confining the system to evolve in a 10 dimensional surface in the full phase space. The 
absence of the BD ghost has been checked from different perspectives and various further 



consistency issues of the dRGT model have also been examined [20M31 



In massive gravity, there is a second fiducial metric, which is promoted to be dynamical 
in bi-metric theories. It has been shown that the bi-metric model generalized from dRGT 
massive gravity is also free of the BD ghost 17, 32|. The multi-Vierbein and multi-metric 
model that are free of the BD ghost have also been formulated (33I . 



34 . Solutions of dRGT 



massive gravity and bi-metric gravity have been investigated by various authors [35H52 



In this paper, we generalize dRGT massive gravity so that the graviton mass can vary 
in time and space. This can be most easily done by introducing a scalar whose background 
value sets the graviton mass. The mass-varying massive gravity model we are proposing 
is free of the BD ghost, as we will explicitly show by examining the constraint system in 
the Hamiltonian formulation. We will also study the FRW cosmology of this model. In this 
scalar-dRGT graviton system, plenty of interesting cosmological scenarios might be possible. 
Among them, we will tentatively discuss a scenario where the scalar plays the role of the 
chaotic inflaton in the early universe and its low energy vev gives rise to dark energy. 

The paper is organized as follows. In Section [Til we describe the model we are proposing. 
In Section IHH we prove that mass- varying massive gravity is free of the BD ghost. We 
will first show that the model has a Hamiltonian constraint and then prove that there is a 
secondary constraint generated by the Hamiltonian constraint and calculate the secondary 
constraint explicitly. In Section ITV] we discuss the homogenous and isotropic evolution of the 
model, focusing on the case when the fiducial metric is Minkowski. In particular, we will de- 
rive the background equations of motion for the case that the dynamical and fiducial metrics 
can be diagonalized simultaneously, and briefly discuss a simple model that unifies inflation 
and dark energy. We conclude and outline possible further consistent generalizations of our 
model in Section IVl 



II. THE MODEL 

As in the Standard Model, the mass of a particle can be generated by a condensate of a 
scalar field. A constant mass can be generated by a vev, while a condensate of the k = 
mode generates a mass parameter that is a function of time. For massive gravity to have 
a varying mass, we may promote the quantity m 2 g M Pl m g being the graviton mass, to be a 
function of a scalar field V(ip). General covariance is spontaneously broken unless the scalar 
sits at V(i/j) = 0. Adding a kinetic term and a potential, we will consider the mass-varying 
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massive gravity model 



S = / d 4 x 



M 2 1 
-^■R + V(iP)(U 2 + a 3 U a + a.U,) - -d^d^ - Wty) 

The dRGT massive gravity 13] is the special case of our model when ip sits in a non-zero 
vev. As in dRGT gravity, a% and are free dimensionless parameters (They may also be 
consistently promoted to functions of ip; see Discussions.), and the three graviton potentials 
read 

U 2 = IC^ U] , V, = KlK?X pV U, = KyC v K^ (2) 
where the antisymmetrization is defined as usual (e.g. /Cf^/C^ = (JC^JC^ — K%Kfy/2) and 



K = % - ^g^f A Bd p <p A dA B - (3) 

Jab is a fiducial metric, which is often chosen as Minkowski Jab = Vab- We will for simplicity 
consider the Minkowski case when we study the cosmology of the model. The four <p A (x) 
are Stiickelberg scalars introduced to restore general covariance (§], so that the theory (pQ) is 
invariant under 

9,u{x) -+ q^q^9 p <t(x), <j) A (x) A (x); x^x'». (4) 

For the case Jab = Vab, the Stiickelberg scalars form Lorentz 4- vectors in the internal space 
(or, field space) labelled by A, B, ... . So in this case the theory also has the (global) Poincare 
symmetry 

<P A ->• 4> A + c A , <P A -> A A B B , when f AB = r) AB . (5) 

We may have a complete gauge fix by conveniently specifying the functional forms of (j> A (x). 
For example, we can choose "unitary gauge" (fi A = x A , leading to /C^ = 5^ — y ' g w f pv . We 
will work in this gauge to show that the theory is free of the BD ghost. Note that, even 
if Jab is not Minkowski, for some cases, we can choose a gauge so that J ? ABd p <fi A d u (j) B is 
Minkowski. 



III. ABSENCE OF THE BOULWARE-DESER GHOST 



In this section, we show that the promotion of the mass parameter to a function of a scalar 
does not spoil the absence of the Boulware-Deser ghost of the dRGT structure. We will first 
show that the Hamiltonian constraint is maintained in mass-varying massive gravity, and 
then prove the existence of a secondary constraint, arising from the time preservation of the 
primary Hamiltonian constraint, and compute the secondary constraint explicitly. We also 
argue that there is no tertiary constraint. To do these, it is more convenient to work in the 



following representation [19 



S 



d 4 xC 



& A X\J—g 



Mi 



R + V(ip) \J2 PneA - \d^d^ - W(iP) 



.n=0 



(6) 
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where 



e = l, e 1 = X» t e 2 = ^^, e 3 = ^;^, e 4 = Xtft v v X p p K\ = det(A^) (7) 
with 

X^y/g^, n = (-l)"Q(4-n)(3-n) + (4-n)a 3 + ^Y (8) 

First, we decompose the two metrics in the ADM formulation to identify the apparent 
dynamical degrees of freedom. Introducing the lapse functions, shift vectors and 3D induced 
metrics, we can rewrite the metrics as 

g^dx^dx" = -N 2 dt 2 + j ij (dx i + N i dt)(dx j + N'dt) (9) 
f^dx^dx" = -L 2 dt 2 + ^ ij (dx i + Ddt)(dxi + tfdt). (10) 

We will lower and raise the spatial indices by %j and its inverse ^ respectively. The fiducial 
metric f^ v , the lapse N and the shift Ni do not have kinetic terms in the Lagrangian (jH]), so 
they are non-dynamical. The potential dynamical degrees of freedom are 7^ and the scalar 
ip. 7jj has 6 components, while a spin-2 field has only 5 propagating modes. If there were 
no constraints on the dynamical degrees of freedom arising from the equations of motion of 
N or N l and their time conservation, then there would be six modes and the sixth mode 
would appear as the Boulware-Deser ghost |8f . dRGT massive gravity, by carefully choosing 
the graviton potential, produces two constraints [If], 17], eliminating the Boulware-Deser 



ghost. In the scalar-dRGT model (0), as we will see, two constraints also arise and the 
Boulware-Deser ghost can also be eliminated. 



A. The Hamiltonian constraint 

Defining the conjugate momenta 

71 1 = , 7T = r, (11) 

and performing the Legendre transformation, we can switch to the Hamiltonian formulation. 
In General Relativity, N and N 1 appear as Lagrange multipliers in the Hamiltonian, so 
their equations of motion clearly impose constraints on the dynamical variables and their 
conjugate momenta. In the model (j6j), none of the ADM variables N and N l appears as 
a Lagrange multiplier, but the 4 equations of motion for N and N l can actually give rise 
to an equation independent of N and N l , i.e., a constraint on 7^, tp and their conjugate 
momenta. That is, there is an implicit Lagrange multiplier in the Hamiltonian. For the case 
at hand, we can make this Lagrange multiplier explicit by re-defining the shift vector N l , 
and the Lagrange multiplier turns out to be N. To achieve this, we introduce a new shift 



vector n 3 , following |10j, [16|, [18 



N l -V = {L8) + ND i j)n j , (12) 
where D l j is defined by the following matrix equation 

xD\D k j = (Y k - D i m n m D k in l )£, kj , x = l-n%r^. (13) 
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The solution to this equation can be written in matrix form as 

D = (xl + nn T £)- 1 y/(xi + nn T ^)r-^, 



(14) 



where the various matrices mean that I = = and T = (jij). With this change of 

variables, after some algebra, the Hamiltonian can be written as 



H = J dx 3 U = J dx 3 (U - NC) 

where 

Ho = —(Ln* + U)TZ t - L^yAV(ip) 
C = TZ + g R i D i j n j + yfyBV(ip). 

Here 1Z, TZi, A and B are given by 

2 ,, I 



f 1 ^ 3 Ml 



IT* 3 ITa - -7T 2 ). 



K=^-(M 2 P ®R - d^&i) - 2W(i/>)) + -^(^ 
2 <fy M J 

A = Pix 1 * + 2 [xD\ + n% J Dhn k ] + (3 3 hx?D* l n%B i k n k + x^D^D 3 ^ + f3 4 
B = (3 + Pix^D\ + i,xi)\i)' , + /3 3 x^D [i i D j j D k] k . 



(15) 

(16) 
(17) 

(18) 
(19) 

(20) 

(21) 



Now, n l are auxiliary fields and can in principle be integrated out by imposing their equations 
of motion dC/dn 1 = —dl-ijdn 1 = 0, or equivalently imposing the following three equations 

C i = U i - V(iP)^x^n% [(3,51 + I3 2 x^{51D\ - D 3 \) + (3 3 x(5jD^ k k D\ - 2D\D>\j 

= 0. (22) 

After this, the equation of motion for N, which is a Lagrange multiplier in the Hamiltonian 
(|T5|) . becomes a constraint 



C{lij, iff, vr, n^Tij, Tr ij , i/>, tt)) = 0. 
This is the Hamiltonian constraint. 



(23) 



B. Existence of a secondary constraint 

Now, a constraint should be preserved under the evolution of the dynamical system, so 
for consistency the following equation should be satisfied 

±C(x) = {C{x), H} = J d 3 y [{C(x), H (y)} - N(y){C(x),C(y)}} = 0, (24) 

where {, } is the Poisson bracket, defined by 



{A(x),B(y)} 



3 , 5A(x) SB(y) 5A(x) 5B(y) 



d*z 



8 lmn {z) 5n mn (z) 5n mn (z) 5 lmn (z) 
f d3z ( 5A(x)5B(y) 5A(x)5B(y) \ 
Stjj(z) 57r(z) 5tc(z) 5i/j(z) J 



(25) 
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If the equation fl2l]) is satisfied trivially, our constraint system is self-consistent and no 
further checking is needed. If not, we should check whether this equation determines the 
lapse N. We can see that if {C(x),C(y)} does not vanish on the primary constraint surface 
C = 0, then the equation (I24p can be used to determine N and no secondary constraint 
arises. If the equation (T23J) can not be used to determine N (and is not trivially satisfied), 
a secondary constraint arises. That is, if {C(x),C(y)} vanishes upon using C = 0, we will 
have a secondary constraint 



C {2 \x) = J d 3 y{C(x),U (y)} = 0. 



(26) 



In the following, we will show that {C(x),C(y)} does vanish upon using C = 0. 

From Eq. ({TBI and n l 's equation of motion, we can infer dC/dn l = 0, which means we 
have 



5C 



DC 



DC 



dir 1 



ii' 



dC 

^^mn _|_ 

dil) 



07T 



5 



7T. 



(27) 



Now, we can expand 

{C(x),C(y)} = {K(x),K(y)} + {KiixlKiiy^D^n^Dhrliy) 



+ 



{n(x),'JZ l (y)}D i k n k (y) + S 



n ' i!! r)^^D i k n k (y)-(x^y) 



57i mn (x) 



where 



S mn (x) 



9lr, 



9lr, 



(28) 



(29) 



In this expansion, we have dropped terms that are proportional to 5 3 (x — y), since 
{C(x),C(y)} is antisymmetrical in exchanging x and y. In particular, the terms involving 
the factor dV(ip)/dip cancel each other. 

The Poisson brackets of TZ and TZi satisfy the same algebra as that of pure General 
Relativity: 



(30) 



{K(x),K(y)} = - [JZ t (x)d x ,5 3 (x - y) - (x y 
{K(x),Ki(y)} = -n{y)d x ,5 3 {x - y), 
{7^- (x), 7^(2/)} = - -R.ji.Di ),. <)'■''(,■ - y )-{xv> y)] 

To show this, one can simply expand the Poisson brackets with 

1 o 



2 



n = n GR - ^(diipd^ + 2W(i/>)) 



2^ 



(31) 
(32) 



utilize the algebra of pure General Relativity [53[, and the smoothing function technique 
introduced below. 

To deal with the Dirac delta function's derivatives, we introduce smoothing functions 
f(x) and g(x), and the quantities 



d 3 xf(x)C(x), G = / d 3 yg(y)C(y). 



(33) 
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We require that f(x) and g{x) fall off quickly enough at infinities so that boundary terms 
can be dropped when integrating by parts. 

So we shall first calculate the Poisson bracket 

{F,G} = J d 3 xd 3 yf(x)g(y){C(x),C(y)}- (34) 

After using the algebra (!30|) . some partial integration and carrying out one of the integrals, 
we get 

{F, G} = - J d 3 x (/ d i9 -gdif)!*, (35) 

where 

V' 1 = (K + 1ZjD j k n k )D\n l + TV + S a - nj D j k n k . (36) 

Then {C(x),C(y)} can be extracted by writing f(x) = J d 3 yf(y)5 3 (x — y) and g(y) = 
J d 3 xg(x)5 3 (x — y) and comparing Eq. (1551) to Eq. (13"4"|) . which gives 

{C(x),C( Z /)} = -[P i (x)^5 3 (x- 2 /)-(x^ 2 /)] . (37) 

To proceed, we have to calculate the partial derivatives d{D\rf) / d^ mn and 8(^8) / d"f mn 
in S mn , using the defining relation of D 1 ^, i.e., Eq. (fT3~j) . For this, we can essentially follow 
the dRGT massive gravity case [16j, as these partial derivations do not involve differentiation 
with respect to ip. After some lengthy calculation, we get 

V i = CD\n l . (38) 

So {C(x),C(y)} vanishes on the surface of C = 0. 



C. The secondary constraint 

Now, we shall derive the expression for the secondary constraint 

C^\x) = J d 3 y{C(x),n (y)}. (39) 

To this end, we can expand the Poisson bracket 

{C(x),U (y)} = -{n(x),n i (y)}(Ln l + V)(y)-D\n\ 

+ VWWA mn (y) + ld\u\x) ^^VWWA mn (y) 

- *~(*) ^ W+l'm - + ^ - v) 

+ L ^ A 8^^ 5 {X ~ V)+ ^ Ifylfr k ( X ~ V) (40) 
where S mn (x) is defined as Eq. ( )29|) and A m ™ is defined as 

A mn = J_ jHvg^) (41) 

V7 ^7mn 
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Note that we have exploited the relations dC/dn l = and dHo/dri 1 = 0. After some 
integration by parts and using the Poisson bracket algebra for 1Z and IZi and Eq. (|34|1 . the 
secondary constraint can be written as 



C (2) = CV l {Ln l + U) + V{4>)L ( lmn n k k - 2vr m „) A mn + 2L^V ' m {V{4>)A mn ) lm D 



+ [RjD\n k - V(^)^ ljk B u ) V l {Ln 3 + V) + (V^ + V i 1l j D s k n k )(Ln i + V) 



8V 

— y^Bdi^Lrt + U) + LAn + L^Ad^D\n k \ 



(42) 



where 



(43) 



D. Absence of a tertiary constraint 



We also have to check whether the time preservation of C^ 2 \ 

^ 2 \x) = {C^(x),H} = 0, (44) 

can give rise to a tertiary constraint. One has to go through a discussion similar to that 
below Eq. f )25|) . But it is easy to see that there is no tertiary constraint. In dRGT gravity, 
the corresponding {C^ 2 \ 1-Lq] and {C^ 2 \C} do not vanish on the constraint surface. As dRGT 
gravity is a limit of mass- varying massive gravity, {C^ 2 \'Hq} and {C^ 2 \C} in mass- varying 
massive gravity also do not vanish on the constraint surface. Therefore the time preservation 
of simply determines the lapse N and no tertiary constraint arises. 



IV. FRW COSMOLOGY 

In this section, we shall study the cosmology of mass-varying massive gravity ([1]) in 
the homogeneous and isotropic background. As mentioned previous, we will choose, for 
simplicity, the fiducial metric to be Minkowski 

Jab = Vab (45) 

for considering the cosmology of the model. Note that this also covers the cases where Jab 
can be brought to the Minkowski metric by general coordinate transformation, as we can 
always choose a gauge for the Stiickelberg fields 4> A . For an interesting cosmology with the 
correct hot history, we shall add some matter source 

S T = S + S m , -^=|%# p!/ = diag(p m ,p m ,p m ,p m ) (46) 
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which have been assumed to be a perfect fluid. For a generic FRW solution, the dynamical 
metric and the fiducial metric do not have to be diagonalized simultaneously (36[. In this 
paper, however, we shall focus on the diagonal case, already with rich cosmologies, leaving 
the non-diagonal case for future work. 

We shall discuss different FRW universes separately. Since our fiducial metric is 
Minkowski, whether we can have simultaneously diagonal metrics for different spatial ge- 
ometries depends on whether the Minkowski metric can be charted by these geometries. The 
flat case can be charted trivially by choosing "unitary gauge" <p A = x A . While the open 



universe can still be charted, the closed universe does not admit a FRW chart 37 



A. Equations of motion 

1. Flat universe 

For this case, the metric and Stiickelberg fields ansatz we are interested in can be written 

as 

ds 2 = -iV(r) 2 dr 2 + a(r) 2 ^dxMx j , 0° = 6(r), ft = a x\ (47) 

where a is constant. Without lost of generality, we may assume that a > 0, a > 0,b' > 
0,N > 0, the action for this background reduces to 



St = I d'.r 

where 



-3M|^ + VW)(v% + a 3 u% + a 4 uf ) + ^' 2 - Na 3 W(^) 



+ S m , (48) 



ul = 3a(a - a )(2Na - a N - ab'), (49) 

uf = (a - a ) 2 (4iVa - a N - 3ab'), (50) 

uf = (a-a ) 3 (N -b'), (51) 

and we have define ' = d/dr . If we also define 



a o . da a 
a Ndr a 



and the polynomials 



h 


l<*i,X) = 


(3 


- 2X) + a 3 (3 • 


-X)(l 


-X) + a 4 (l 


h 


>j,X) = 


(1 


-X) + a 3 (l- 


xf + 




h 


[an,X) = 


(3 


-X) + a 3 {l- 


X), 





X)\ (53) 
(54) 
(55) 

then the equations of motion for 6, N, a and if) are given by, respectively, 

VWHMan, X)+V(i))f 2 (a i , X) = 0, (56) 
3M 2 P H 2 = p MG + Pm , (57) 

-2MpH = puG + PMG + Pm+ Pm, (58) 

dW dV 

^ + 3Hilj + — + —i(X- l)Cf 3 K X) + h{a h X)) + 3&f 2 K X)] = 0, (59) 
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where p m and p m are the density and pressure for the matter source respectively and the 
effective density and pressure are given by 

Pug = + W(i/>) + Vty)(X - l)h{on, X) + V($){X - X), (60) 

Pmg = \i> 2 - W(if>) - VW)(X - l)h{a h X) - VW){b - l)/iKX). (61) 

Note that /i(aj,X) and /^(a^X) satisfy a differential relation 

x] d[3X- 3 / 2 (^,X)] d^M) 
/i = = 3f 2 (a i7 X)- . (62) 

Eq. ( 15 6 p can be integrated by the method of viable separation, which reduces to an algebraic 
equation for V{jp) and a, 

V(^) = 3V e~ J ^ da = 3 , a ° y ° = const. > 0. (63) 

Gr7 2 (ai,X) 

That is, this equation can be seen as a constraint. 

Unlike the diagonal case of dRGT massive gravity, the system (I56H59I) can have interesting 
expanding solutions. But the system is highly non-linear, and, apart from some simple ones, 
exact solutions are difficult to find. One good strategy might be to apply a dynamical system 
analysis on it to extract its phase space information (see e.g. For this, one shall first 

reduce it to a standard dynamical system and analyse the equilibria the system possesses, 
which we leave for future work. 



2. Open universe 



For an open FRW universe, K < (we shall use the constant k = y/\K\), the metric and 
Stiickelberg fields ansatz we are interested in can be written as 



ds 2 = -X(r) 2 dr 2 + a(r) 2 %dxW, fiydaW = S ij dx i dx j 



l J rk 2 (5 ij x i x^ l 



<f = b(r)Jl + WiSijxW), <p l = kb(r)x\ 



(64) 
(65) 



For this ansatz, g^ u and rjABd^^dycj) 3 become diagonal simultaneously, so )C£ becomes 
diagonal and its square root is uniquely defined and easy to calculate. Assume that a > 
0, b > 0, b' > 0, X > 0, the action now reduces to 



Si 



- 3M l P (k 2 Na-\ 



a' 2 a. 



o 



.On 



X 



a 



+ - Na 3 W&) 



where 



u 



o 



u 



o 



u 



o 



3a(a - kb)(2Na - kNb - ah') 
(a - kbf{ANa - kNb - 3ab') 
(a - kbf(N - b') 



(66) 

(67) 
(68) 
(69) 
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Defining 

. da a ,„ n . 

Y=—, a=—-, H=~ 70 
a IS (xt a 

the equations of motion for b, N, a and if) are given by, respectively, 

Vty) ( H -^j /i(a^) + ^)/2(«„y) = 0, (71) 
3fc 2 

M|(3tf 2 -~-j-) = Pmg + An, (72) 
k 2 

M 2 P {-2H - 2—) =p MG + p MG + p m + p m , (73) 
dW dV 

i/j + 3Hi; + — + —[(Y- l)(/ 3 (ai, V) + A (a*, F)) + 36/ 2 (a l; F)] = 0, (74) 

where the polynomials fi{oLi, Y), /^(a:*, F) and fs(cti, F) are defined the same as that in the 
flat universe case, p m and p m are the density and pressure for the matter source respectively, 
and the effective density and pressure are given by 

Pmg = \i? + Wty) + V(^)(Y - 1)/ 3 K F) + VtyW - l)f x {a h F), (75) 

Pmg = ~ W(i/>) - VU>)(Y ~ 1)/ 3 K, Y) - V{^)(b - l)h{a h Y). (76) 

Again, in general, the system (I71ti74p is highly non-linear, and, exact cosmological solutions 
are difficult to obtain and one may apply a dynamical system analysis on it to extract its 
phase space information. 



B. A special solution and its cosmological implications 

In the open universe case, there exists one simple branch of solutions if we have 

a A = — , Y = — = , (77) 

4 a a 3 

that is, if we impose a constraint on the free parameters a 3 and a A and let the two scale 
factors b(t) and a(t) scale with each other. This gives rise to 

h(a i ,Y) = f 2 (a i ,Y) = 0. (78) 

In this branch, the dynamical system (171ti74p reduces to that of a General Relativity plus a 
scalar 

3fc 2 1 

Mp(3H 2 --^) = + WW) + Pm , (79) 
k 2 

M 2 (-2H-2—)=i; 2 + p m +p m , (80) 

(X 

$ + 3HtP + ^j- = 0. (81) 
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where 

W{$) = W{$) - \vty). (82) 

The evolution of this system will determine the effective mass ^/V(ip)/Mp of the graviton. 

As an aside, we could have for simplicity set W(ip) to so that the scalar potential 
W(ip) would be simply supported by V(ip). For the graviton sector to be free of tachyonic 
instabilities, we need V(ijj) > 0, which leads to W(ip) < 0, that is, the scalar system can 
only have negative potential energy (we can still require it to be bounded below though.). 
In this case, the graviton mass could only increase and we would not have a de Sitter-like 
solution but have an anti-de Sitter-like solution instead. 

In this special case, the background evolution of mass-varying massive gravity connects 
to some well studied cosmological systems. We may identity this scalar with inflaton, 
quintessence, quintessential-inflaton, or any other scalar that plays a role in the homo- 
geneous and isotropic evolution of the universe. In constructing these scenarios, however, 
one shall make sure the effective graviton mass in the current cosmic time is small enough 
to satisfy all the existing General Relativity tests. 

For example, here we discuss a simple model that gives rise to inflation and dark energy. 
We shall consider 

W(i/>) = M 2 ifj 2 , V(i/>) = ^a 2 3 mM 2 {ifj - m), (83) 

which leads to an effective potential W = M 2 (ip — m) 2 + m 2 M 2 , where m and M are free 
parameters of the model, to be adjusted to satisfy observations. To give a concrete example, 
if we choose M ~ 10~ 5 M P and m ~ W 5 H (H is the current Hubble constant), and set the 
ip field initially to be around the Planck scale Mp, then the universe will start with chaotic 
inflation, reheat around the minimum of W at ip = m ~ 10 5 H and has, at the background 
level, an effective cosmological constant m 2 M 2 ~ H^Mp, which gives rise to dark energy. 
Interestingly, around this vev, the graviton is massless and its potential starts only from the 
cubic order in perturbation. One may further consider whether this model will reproduce all 
the successes of inflation's power spectra. For this one shall study the perturbations around 
the quasi-de Sitter background generated by W(ip). Here we give a simple argument which 
implies that these successes might be reproduced. In massive gravity, there is a Vainshtein 
radius for some matter source, within which massive gravity effectively reduces to Einstein 
gravity. In the early universe the squared mass of the graviton is ~ a 2 mM 2 ip/Mp. For 
a spherical source M. the Vainshtein radius for mass-varying massive gravity is given by 
(M./m 2 Mpy/ 3 , where m g is the graviton mass. We may see a Hubble patch as a spherical 
source M. ~ pxH~ 3 ~ MpH 2 ■ H~ 3 (H' 1 is the Hubble radius) and estimate the Vainshtein 
radius as 

r v = (a 2 3 H H 2 r l *. (84) 

Even if the free parameter a 3 is chosen as 0(1), the Vainshtein radius r v is still much 
greater than the Hubble radius at that time. Therefore we may expect that all the success 
of the conventional inflation model will still hold. However, this should be confirmed by an 
explicitly perturbation analysis, which might also give rise to some testable differentiating 
effects. We leave this for future work. After reheating the inflaton settles down on the vev of 
the effective potential, which is of the dark energy scale, and the graviton becomes massless, 
so the big bang hot history is recovered afterwards. 
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V. DISCUSSIONS 



Particle masses can often be generated by a scalar field, like Higgs. In this paper, we have 
introduced a model where the graviton mass is effectively generated by a scalar field, which 
spontaneously breaks general covariance of General Relativity. This model is a generalization 
of dRGT massive gravity and is also free of the BD ghost. We have proven the absence of the 
BD ghost by explicitly obtaining two constraints in the Hamiltonian analysis of the theory. 

We may further generalize our model by promoting 0:3 and 04 to functions of if), i.e., 
otz{ip) and a^ip). It is straightforward to check that a similar proof can go through, thus 
the BD ghost is not re-introduced either. We may also consider a non-minimal coupling 
y/—gh(i())R, h(ip) being a general function of ip. For this case, we can do a conformal 
transformation, which reduces the non-minimal coupling to the Einstein case but introduces 
some non-canonical kinetic terms and scale the graviton potential terms with powers of h{ip). 
We also expect this generalization is free of the BD ghost. Further more, we can make the 
fiducial metric a genuine dynamical one 17| . which should also not spoil the absence of the 
BD ghost. 

We have also initiated the study of the homogenous and isotropic cosmology for this 
model. We have derived the background equations of motion for the case when the two 
metrics can be diagonalized simultaneously, pending a detailed analysis of the systems of 
differential equations and the non-diagonal case for future work. Nevertheless, we have 
found a simple solution when and a± satisfy the relation f J77|) . In this special case, the 
model's background evolution completely reduces to that of General Relativity plus a scalar. 
We have tentatively considered a simple scenario where the mass- varying scalar is identified 
with the chaotic inflaton and its low energy vev gives rise to the current cosmic acceleration. 

In that case, the graviton mass is initially relatively large and settles to zero at reheating. 
Of course, one may alternatively consider various other scalar potentials, many of which 
might also lead to interesting cosmologies. If we draw more direct comparison with the 
spontaneous symmetry breaking in the Standard Model, one may consider a scalar potential 
where the graviton is initially massless and gains a small mass after symmetry breaking. 
Then, an interesting question may be: what is the relation between this scalar and the 
Higgs boson (or the Higgs bosons)? A scalar coupled to gravity might respond substantially 
to some local matter distribution, which might induce spatial variation of the graviton mass. 
As we generally require the current graviton mass is small, the spatial variation effect may 
arguably be small, but it is interesting to understand whether the current experiments and 
observations can impose any theoretical bounds on mass- varying massive gravity from this 
perspective. 
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